Entanglement dynamics of three-qubit states in noisy channels 
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We study entanglement dynamics of the three-qubit system which is initially prepared in pure 
Greenberger-Horne-Zeilinger (GHZ) or W state and transmitted through one of the Pauli channels 
<Tz, o~ x , o y or the depolarizing channel. With the help of the lower bound for three-qubit concurrence 
we show that the W state preserves more entanglement than the GHZ state in transmission through 
the Pauli channel a z . For the Pauli channels a x , a y and the depolarizing channel, however, the 
entanglement of the GHZ state is more resistant against decoherence than the W-type entanglement. 
We also briefly discuss how the accuracy of the lower bound approximation depends on the rank of 
the density matrix under consideration. 

PACS numbers: 03.67.Mn, 03.67.Hk, 03.65.Yz. 
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I. INTRODUCTION 

Entanglement is the fundamental property of quan- 
tum systems which has been found as a resource for 
many potential applications, such as quantum tclepor- 
tation, cryptography and superdense coding [l| , to name 
just a few. Besides the other manifestations of entan- 
glement, quantum teleportation still holds central place 
in quantum information theory, due to a few remarkable 
recent experimental achievements 0-0] • A protocol for 
quantum teleportation was originally suggested in the 
seminal works of Ekert [5[ and of Bennett et al. [6[ to 
facilitate communication between two partners. The pro- 
tocol enables one to transmit an unknown quantum state 
to the remote recipient by means of an initially shared 
two-qubit entangled (Bell's) state. Soon thereafter, com- 
munication protocols for quantum teleportation between 
several partners, that are based on many-qubit entangled 
states, were suggested @,Hj]. In particular, Karlsson and 
Bourennane [7| suggest protocols for quantum telepor- 
tation between two and three partners with three-qubit 
GHZ state 
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(|000) + |111)) 



(1) 



Later Agrawal and Pati suggested using three-qubit W 
state 



|W0 = 1(V2|001) 



|010) + 1 100) 



(2) 



in protocols for quantum teleportation between two part- 
ners and superdense coding |8J. Nowadays, protocols for 
quantum teleportation, that are based on many-qubit 
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entangled states, are widely used in order to construct 
quantum networks [9l-[Tlj|. 

While protocols for quantum teleportation require 
pure entangled states, in practice one has to deal with 
entangled states that are mixed due to imperfect control 
during their transmission through communication chan- 
nels. It is known that a pure entangled state can be ex- 
tracted from several copies of mixed entangled states, if a 
so-called protocol for en tang lement purification is applied 
to the mixed states pI2Hl4l j . To provide the protocol on 
received mixed states it is necessary to know how much 
entanglement is preserved in a single mixed state after 
transmission through a (noisy) communication channel. 

The entanglement dynamics of many-qubit entangled 
states under influence of different types of environment 
has been intensively discussed during last decade [l5rll8l | . 
For instance, entanglement dynamics of three-qubit GHZ 
and W states for special types of environment coupling, 
such as a thermal bath at zero and infinite temperatures 
as well as the dephasing (i.e. the Pauli a z ) channel, was 
discussed in Ref. [l5j]. In particular, with the help of the 
lower bound for three-qubit concurrence it was shown 
that GHZ-typc entanglement is more fragile under these 
types of environment coupling than the entanglement of 
the W state. 

In this contribution we also study the entanglement dy- 
namics of three-qubit GHZ (JXJ) and W © states under 
the influence of its environment. We compliment results 
which were obtained in Ref. fl5j | by considering entan- 
glement evolution of the states ([Ill-Q in transmission 
through various noisy channels, namely, the Pauli chan- 
nels a z , a x and a y as well as the depolarizing channel 
Although the entanglement dynamics of the three-qubit 
states in the Pauli channel a z was considered earlier in 
Ref. [l5j we include this channel in our discussion in order 
to complete the investigation of the Pauli channels. 

The time evolution of three-qubit (open) quantum sys- 
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terns which are initially prepared in the pure states (Q}- 
@ and transmitted through one of the Pauli channels or 
the depolarizing channel was recently analyzed by Jung 
et at (l9|. For each of these channels, the time evolution 
of the GHZ and W states is given by mixed state den- 
sity matrices which were obtained by analytical solving of 
corresponding master equations in Lindblad form. With 
the help of the analytic expressions for the mixed state 
density matrices we describe the entanglement evolution 
of the mixed GHZ and W states. To quantify the entan- 
glement of the mixed states we use the lower bound for 
the three-qubit concurrence which was recently suggested 
by Li et aZ.[2(|. We show that the W state preserves 
more entanglement than the GHZ state in transmission 
through the Pauli channel a z . Surprisingly, for the Pauli 
channels o~ x and o~ y as well as the depolarizing channel 
GHZ-type entanglement is more resistant against deco- 
herence than W-type entanglement. 

In addition, we briefly discuss how accuracy of the 
lower bound approximation for three-qubit concurrence 
depends on the rank of the density matrix. In general, 
to quantify entanglement of a three-qubit mixed state, 
the convex roof extension for three-qubit concurrence is 
already known 21]. Unfortunately the calculation of the 
convex roof implies an optimization procedure which still 
has no analytic solution. Of course, the optimization 
can be done numerically. However, the calculation of 
the convex roof for a density matrix with rank r leads 
to optimization over r 3 free parameters [l6j which is in- 
deed quite a formidable task. In practice moreover, one 
is mostly interested in the minimal (nontrivial) amount 
of entanglement which is preserved in a mixed state. 
With regard to this practical requirement several ana- 
lytical lower bound approximations for the concurrence 
were recently suggested [l(| IHJ and were justified 
to be proper entanglement measures. It is important to 
know how accurate the approximation is and from which 
parameters of the density matrix this accuracy depends 
on? As we have already mentioned in this work we use 
the lower bound approximation for the concurrence from 
Ref. [20( in order to describe entanglement evolution of 
the mixed states (pTJ)- (ED - We found that for the density 
matrices with rank r < 4 the lower bound coincides with 
the convex roof for the concurrence. For the density ma- 
trices with higher rank the application of the lower bound 
becomes restricted: the lower bound can not describe 
long-time evolution of the mixed state. 

The paper is organized as follows. In the next sec- 
tion we present the necessary theoretical background to 
describe the time evolution of the states (HJ-© in the 
Pauli channels o~ z , o~ x and o~ y as well as the depolarizing 
channel. In Subsection III Bl we discuss the convex roof 
extension for many-qubit concurrence and introduce the 
lower bound to be used. With the help of presented tech- 
nique we show the entanglement dynamics of the GHZ 
and the W states in transmission through the Pauli and 
the depolarizing channels in Sections Mil and IIVI respec- 
tively. In Section [V] we summarize the results which were 



obtained in the previous two sections. We also briefly dis- 
cuss the relation between the accuracy of the lower bound 
approximation and the rank of the density matrix. 

II. THEORY 

A. Time evolution of states transmitted through 
noisy channels 

Decoherence is a well-known quantum phenomenon 
that occurs for all open systems, i.e. if they are cou- 
pled to their environment. In order to classify and better 
understand this (often undesired) interaction of a system 
with its surroundings, a large number of noise models 
has been investigated during the last decades [H, [24| , in- 
cluding various thermal bathes at either zero or infinite 
temperature, the phase and the amplitude damping, or 
just the Pauli channels <Jn_ (a = x,y,z), etc. Following 
similar lines, Jung et al. [19( have recently analyzed the 
time evolution of the three-qubit GHZ (TJJ and W © 
states, if they are transmitted through one of the Pauli 
or the depolarizing channel, in order to explore the effi- 
ciency of the 'two-sided' teleportation protocols that are 
based on these entangled states. In this work [l9| . in 
more detail, an initially pure entangled state p(0) was 
supposed to be transmitted through (one of) these chan- 
nels for the time t, and its time evolution p{t) obtained 
as solution of a (Lindblad-type) master equation 

^ = -i[H s ,p]+J2 (u,apLl a ~ \{L\ M U^p\\ .(3) 

i,a 

In this master equation, the (Lindblad) operators Li. a 
were assumed to act independently upon the i— th qubit; 
for example, the operator L\ iZ = \fko z ® 1 ® 1 describes 
the decoherence of the first qubit under a phase-flip er 2 , 
and where the coupling constant k is approximately in- 
verse to the decoherence time with regard to such a phase- 
flip. Later we shall also refer the Pauli channels o~ x and o~ y 
to bit-flip and bit-phase-flip coupling of the three-qubit 
system to the environment |l[ . For any given Pauli chan- 
nel o~ a , therefore, the master equation Q only includes 
three Lindblad operators, Li, a , L2, a and L^, ai while nine 
of these operators are needed for the depolarizing chan- 
nel, Li t0l , (i = 1,2, 3, a = x,y,z). In the latter case, 
each of the qubits can be affected with equal coupling 
strength by all three Pauli channels simultaneously [l|. 

B. Concurrence for three-qubit mixed states 

Knowing the time evolution of the three-qubit GHZ 
(fT]) and W (J5J) states in transmission throw the Pauli and 
the depolarizing channels, we still need to quantify the 
remaining entanglement of the mixed states after they 
have been passed through some of these noisy channels 
discussed above. For mixed states, in fact, any exact 
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quantification of their entanglement has been found dif- 
ficult, and no general solution is known until now [2l[ 
apart from two-qubit systems. In the latter case, Woot- 
ters concurrence [25| provides a very powerful measure 
but, despite of its rather simple form for just two qubits, 
there is not unique generalization of this measure avail- 
able even for mixed bipartite states, if the dimensions of 
the associated Hilbert (sub-) spaces are larger than two 
[22l |H, [27j ■ The generalization of the concurrence for 
mixed many-partite states still remains an open problem. 
However, several useful lower bounds for the concurrence 
have been proposed in the literature [H, EE Hfll HI] m 
order to characterize the entanglement of many-partite 
states. Here, we shall not discuss these suggested 'mea- 
sures' in further detail but make use of a recent work 
by Li et al. [20| who suggested the lower bound for the 
concurrence for three-qubit states which is based on the 
positive partial transpose and realignment separability 
criteria. 

Before we shall further discuss the entanglement of 
mixed states, let us first consider the concurrence of some 
pure state \tp). For any (pure) three-qubit state, for ex- 
ample, the concurrence can be expressed as [20] 



CM) = sj\ (3 - Trpf - Trp 2 - Trp§) , (4) 

and where the reduced density matrices pi = Tr 3 fc \ip) (ip\ 
with i =/= j k are obtained by tracing out the re- 
maining two qubits. Using the definition Q, we can 
easily calculate the concurrence C 3 (|GHZ)) = and 
C 3 (|W)) = v/378 for pure GHZ © and pure W © states 
which are considered in the present work. Although both 
states are known to be fully entangled [2lj], t ne defini- 
tion (U) results in two different values; therefore, we shall 
re-normalize the expression (|4]) for each state in such a 
way, that we have C* 3 (|GHZ)) = C 3 (|W)) = 1. This re- 
normalization is justified from experimental viewpoint, 
since a three-qubit maximally entangled state can be 
viewed as a single unit of a quantum teleportation pro- 
tocol. 

Of course, any mixed state can be expressed also 
as a convex sum of some pure states {l^j)} 1 P — 
J2i Pi (V'il- Using this representation, definition (|4]) 
for the entanglement of pure states can be generalized for 
mixed states and gives rise to the so-called convex roof 
(extension) [2l| 



Cs(p) = mm}^pi C 3 (\ipj)) 



(5) 



where the minimum has to be found with regard to all 
possible decompositions of p into pure states \tpi) (and 
with positive coefficients pt). The maximal number z max 
of pure states (in the ensemble) is called the cardinality 
of the decomposition. This cardinality is not fixed by the 
rank r of the density matrix, though it is usually assumed 
to be i max < r 2 [l||[28|]. In practice, therefore, one would 
need to search for all decompositions of p into up to r 2 



pure states in order to just evaluate the convex roof for 
a three-qubit density matrix, implying an optimization 
procedure with r x r 2 ~ r 3 free parameters [la ]. This is 
indeed quite a formidable task already for three qubits, 
and the situation becomes worse since no numerical al- 
gorithm could guarantee to find the global minimum for 
the expression on the rhs of Eq. ([5]). An analytical so- 
lution is known for this optimization task but refer to a 
very special case [29| . 

Instead of finding the exact minimum for the three- 
qubit concurrence ((5]), the rather simple lower bound to 
this measure was suggested by Li et al. [20] 



\ 
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(cf 12 ) 2 



(C 23|1 )2 



(6) 



i=l 



which is given in terms of three bipartite concurrences 
that correspond to possible bipartite cuts of the three- 
qubit system. The bipartite concurrence was originally 
introduced in Ref. [30j and was proved to be an entan- 
glement measure. The bipartite concurrence C 12 ' 3 for 

12 1 3 

qubits 12 and 3 is given by sum of six terms C i in 
Eqn. ((S]); each term is expressed as 



Cf |3 = max{0,Ai-A 2 -A 3 -A4}. 



(7) 



where the A^, k = 1..4 (for given i) are the square roots 
of the four nonvanishing eigenvalues in decreasing or- 
der of the matrix p p] . These matrices are not her- 
mitian and are formed by the density matrix p and its 
complex conjugate p*, and which is further transformed 

by the operators {S^ 3 — Lj 2 ® Lq, i = 1...6} as: 

p* 2 ' 3 = S^ 12 ' 3 p* iS* 2 ' 3 . Details about the construction 
of these operators and the mathematics behind can be 
found in Ref. 20]. In this notation, moreover, Lq is 
the single generator of the group SO (2), while the Lj 2 
are the six generators of SO (4). When the generator Lq 
simply equal to the second Pauli matrix cr yi the gen- 
erators Lj 2 can be expressed by means of the totally 
antisymmetric Levi-Cevita symbol in four dimensions, 
i.e. (L k i) mn = -iskimn, k,l,m,n = 1..4 [HI, HI]. The 
bipartite concurrences C 13 ' 2 and C 23 ' 1 are defined in the 
same way as above. 

With these remarks about the lower bound for the 
concurrence we now have two possibilities to de- 
scribe the time evolution of the entanglement C 3 (p(i)) 
if an initial state p(0) is passed through a noisy chan- 
nel. To determine the convex roof ((5]) as function of 
time t, we would first need to split time t into a set of 
steps, < ti < ... < tiv, and optimize at each time step 
Cs{p{ti)) with regard to ~ r 3 parameters, where r is the 
rank of the density matrix. Finally, an interpolation for 
all times 'in between' these steps need to be performed. 
If, in contrast, we make use of the lower bound ([6]), we 
can evaluate this bound for the density matrix p(t) ana- 
lytically for all times t. In the next two sections, we shall 
therefore apply the lower bound ^ to the concurrence 
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Ts(p(t)) to analyze the decay of the entanglement for an 
initially pure GHZ ((T|) and pure W state © separately. 
However, since the lower bound ([6]) is only an approxima- 
tion for the convex roof ([5]) , we shall explore the validity 
of this approach if the density matrix p(t) departs more 
and more from a pure state due to its interaction with 
the noisy channels. In rather simple case when the den- 
sity matrix has rank r < 4 we shall compare the lower 
bound for the concurrence with actual value of the con- 
vex roof obtained numerically. In Section [V] we discuss 
how the accuracy of this lower-bound approximation is 
related to the rank of the density matrix which differs in 
dependence of the initial state and considered channel. 



III. ENTANGLEMENT EVOLUTION OF AN 
INITIAL GHZ STATE UNDER NOISE 

A. Pauli channel a z 

If an initially pure GHZ state ([I} is transmitted 
through the Pauli channel a z , its time evolution is ob- 
tained as solution of the master equation (J3J with Lind- 
blad operators (Li >z , L% >z , La jZ ) and can be expressed in 
terms of the rank-2 density matrix 19] 



with 



Pit) = \ (1000) (000| + |111) (111|) 



(8) 



-6kt 



(|000) (111| + |111) (000|) 



For this mixed state, the lower bound ^ to the three- 
qubit concurrence is a monoexponcntial function of time, 



t 3 (p(t)) = e 



-6kt 



(9) 



Since the rank of the density matrix ([8]) is two, the convex 
roof extension ([5]) for this density matrix can be even cal- 
culated analytically [15j . In this case, the convex roof is 
shown to follow the behavior of the nondiagonal elements 
(up to the normalization factor). In fact, the convex roof 
for the density matrix ((5J) coincides with the lower bound 
©• 



B. 



Pauli channel <r. 



If the GHZ state (JIJ is instead transmitted through 
the Pauli channel <j x , its time evolution is given by the 
rank-4 density matrix [l9T | 



P(t) 



( a+ 




















a+ \ 





a_ 














a_ 











a_ 








OL- 

















a_ 


a_ 




















a_ 


a_ 

















a_ 








OL- 











o_ 














a ■_ 





V a + 




















a + J 



(10) 



l + 3e 



-4kt 



and 
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In this case, the lower bound (j 
currence becomes 

ts (p(t)) = e 



to the three-qubit con- 



-4kt 



(11) 



i.e. the entanglement of the initial state decays less 
quickly for a bit-flip coupling of the three qubits to the 
environment than for a phase-flip. 

For the rank-4 density matrix (fT0|) we also calculated 
numerically the convex roof ([5]) . We found that the lower 
bound (fTTj) coincides with the numerical values of the 
convex roof. 



C. Pauli channel a y 

For a transmission of the GHZ (Q} state through the 
Pauli channel a yi the density matrix [19J 
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Pit) = 3 



has full rank (i.e. rank 8 

01 = 3e- 2kt + e- 6kt 



with the two functions 



(12) 



and /?2 = 



~2kt 



-6kt 



respectively. For this matrix, the lower bound ([6]) to the 
concurrence gives rise to 



t-a (P(t)) = max{0, i ( 



o — 2kt , —Akt , 

Se + e + e 



l)}, (13) 



or, in other words, this lower bound vanishes already 
after some finite time. Us ing the positive partial trans- 
pose separability criteria 21], we verified however that 
the state (|T2|) becomes separable only asymptotically for 
t — y co , which implies that the lower bound (|13[) does 
not describe the long-term behavior of the entanglement 
of an initial GHZ state if its is affected by bit-phase-flip 
noise. 

For the rank-8 density matrix (|T2l the numerical cal- 
culation of the convex roof ([5]) requires optimization over 
8 3 = 512 free parameters. The numerical value of the 
convex roof ([5]) for the rank-8 density matrix (|T2|) as well 
as for other rank-8 density matrices discussed below has 
not been obtained by us yet. 



D. Depolarizing channel 

If the state (JTJ) is transmitted through the depolarizing 
channel, its density matrix has also rank-8 and takes the 
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T3(P) 




0.2 0.4 0.6 0.8 1.0 

FIG. 1: (Color online) The lower bound (|6} for the three- 
qubit concurrence t 3 as function of time t for an initial GHZ 
state |T]), if transmitted through various noisy channels: Pauli 
channels a z (solid red), a x (dashed green), a y (dotted blue) 
and the depolarizing channel (solid black). 



to the Pauli channels a z . The depolarizing coupling of 
the three-qubit system to the channel is the most de- 
structive for the entanglement. It is also remarkable that 
for density matrices with rank-2 and rank-4, the lower 
bound coincides with the convex roof and describes the 
entanglement evolution for all times, while this bound 
is not applicable for the long-time description of density 
matrices with rank-8 (the Pauli a y and the depolarizing 
channels) for which it vanishes at a finite time. 
IV. ENTANGLEMENT EVOLUTION OF AN 
INITIAL W STATE UNDER NOISE 

A. Pauli channel a z 

A similar analysis as in Section Hill can be made if the 
system is initially prepared in a W state. If the state ^ 
is transmitted through the channel a z , its time evolution 
is described by the rank-three density matrix [19J 
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and this gives rise to the lower bound 
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with 
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l + 3e" 



-Skt j 

e and 7 



4e" 



Here, again, the lower bound ^ to the entanglement 
vanishes already after some finite time due to the condi- 
tion 



r, </>!/» max-'i), 1 (4e~ 12fct + e~ 8kt - l)} . 



(15) 



Fig- HI displays the time-dependent lower bound ((6|) for 
initial GHZ state ([1} if it's transmitted through the dif- 
ferent channels. In all cases, this lower bound decays 
exponentially due to the noise of the channel; in trans- 
mission throw the Pauli channels a x and a y the entan- 
glement of the GHZ state decreases slowly comparing 



-4kt 



(17) 



for the evolution of the entanglement, which moreover 
coincides with the convex roof ([5]) as we verified numeri- 
cally. 



B. 



Pauli channels a x and a v 



If the (initially prepared) W state is transmitted 
through the Pauli channels a x or a v , a full rank-8 density 
matrix is obtained for its time evolution 19] 
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and where the + sign refers to the a x and — to the a y 
channel, respectively. The time-dependent parameters in 
expression (|18|) are given by 
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Since two density matrices p(t)± have the same structure 
of matrix elements, the lower bounds for these density 
matrices coincide. Unfortunately, achieved analytic ex- 
pression for the lower bound for the density matrix (|T5|) 
has no compact form and, thus, we do not show it here 
explicitly. At Fig. [5] the lower bound is shown with blue 
dashed line. As for all rank-8 density matrices above the 
lower bound for the density matrix (|18|) vanishes after 
finite time. 



T3 (P) 




FIG. 2: (Color online) The same as in Fig. [TJbut for an initial 
W state the lower bound © for the three-qubit concur- 
rence T3 as function of time t is shown for noisy channels: 
Pauli channels a z (solid red), a x and a y (dashed blue) as well 
as the depolarizing channel (solid black) 



C. Depolarizing channel 

Finally, if the W state @ is transmitted through the 
depolarizing channel, the density matrix p{t) has also 
rank-8 and is given by 19] 
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where 
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The time-dependent lower bound ((6|) for initial W state 
(|2|) transmitted through the different channels is shown at 
Fig. [5J As in the case of the GHZ state the lower bounds 
for the W state decay exponentially due to the noise of 
the channels. In contrast to the GHZ state, the entan- 
glement of the W state decreases slowly in transmission 
throw the Pauli channel a z comparing to the Pauli chan- 
nels a x and a y . However, the depolarizing coupling of 
the three-qubit system to the channel is again the most 
destructive for the entanglement. For the rank-3 den- 
sity matrix, moreover, the lower bound coincides with 
the convex roof and describes the time evolution of the 
entanglement for all times, while this bound is not suit- 
able for the long-time description of density matrices with 



rank eight (the Pauli a x and a y as well as depolarizing 
channels). 



V. RESULTS AND DISCUSSION 

In previous two sections we showed explicitly the en- 
tanglement evolution of the three-qubit system which is 
prepared in pure GHZ ([T]) or W ([2]) state and transmit- 
ted for the time t through one of the Pauli channels <7 Z , 
a Xl a y or the depolarizing channel. Having these results 
we can investigate entanglement of which state between 
the GHZ and the W is more resistant against decoher- 
ence in transmission through the noisy channels. For the 
Pauli channel a z the lower bounds for the GHZ and the 
W states are given by Eqns. (j9]) and (fT7|) respectively 
and are shown with red solid lines at Figs. Q][2] In fact, 
the W state preserves more entanglement than the GHZ 
state in transmission through the Pauli channel a z for all 
times t. This result was obtained earlier in Ref. [l5|]. For 
the Pauli channels a x , a y and the depolarizing channel, 
in contrast, the lower bounds for the GHZ state, that are 
shown at Fig. [I] are always higher than corresponding 
lower bounds for the W state at Fig. [2j The entangle- 
ment of the GHZ state is thus more resistant against 
decoherence than W-type entanglement in transmission 
through the Pauli channels a x , a y and the depolarizing 
channel. 

Our result extend the investigation of the entangle- 
ment dynamics of the three-qubit GHZ and W states in 
transmission through noisy channels, which was started 
in Ref. [15j- As it has been already known, the W state 
preserves more entanglement than the GHZ state if it is 
coupled to the thermal bath at zero or infinite tempera- 
ture, or the Pauli channel a z [HI]. We showed that the 
W state preserves less entanglement than the GHZ state 
if it is coupled to the Pauli channels a x or a v , or the 



7 



depolarizing channel. 

To describe entanglement evolution of the mixed GHZ 
and W states we used the lower bound ([6]) which is an 
approximation for the convex roof for three-qubit con- 
currence ([5]). We showed that the lower bound coincides 
with the convex roof for the density matrices with rank 
r < 4. This statement, however, is true only for the den- 
sity matrices which were considered in Subsections IHI Al 
IIIIBI and llV Al Whether the lower bound coincides with 
the convex roof for an arbitrary density matrix with rank 
r < 4 introduces an open question of great importance. 
We like to investigate this question in the nearest future. 

For the rank-8 density matrices discussed in this work 
the lower bound for the concurrence vanishes after the 
finite time. Vanishing of the lower bound, however, does 
not mean vanishing of the entanglement of the mixed 
states. We verified with the positive partial transpose 
separability criteria (2l[ that the mixed states become 
separable only asymptotically for the limit t — > oo. For 
the rank-8 density matrices the lower bound for the 
concurrence can describe entanglement evolution of the 
mixed state for a restricted time t < to only, where to 
denotes the time when the lower bound vanishes. The 
comparison of the lower bound with the convex roof for 
the rank-8 density matrices has not been yet performed, 
since the calculation of the convex roof introduce quite 
difficult problem which requires the numerical optimiza- 



tion over 8 3 = 512 free parameters. This comparison is 
left for future research. 

The lower bound approximation for the concurrence 
introduces a powerful simple tool to describe the entan- 
glement of an arbitrary N-partite mixed state. System- 
atic investigation of the accuracy of the approximation 
has not already been done. The fact that the accuracy 
of the lower bound depends on the rank of the density 
matrix is intuitively understandable, since the complex- 
ity of the calculation of the convex roof also depends on 
the rank. In this work we started investigation of the ac- 
curacy of the lower bound approximation with regard to 
the rank of the density matrix. On particular examples 
we showed that the lower bound ©: i. coincides with 
the convex roof © for the density matrices with rank 
r < 4; and ii. vanishes after the finite time for the rank-8 
density matrices. In future we like to investigate in more 
details how the accuracy of the lower bound approxima- 
tion for three-qubit concurrence depends on the rank of 
the density matrix. 
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